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1.0  INTRODUCTION 


A  considerable  amount  of  theoretical  and  numerical  analysis  has  been 
published  recently  in  the  general  area  of  unsteady  flows  and,  in 
particular,  on  unsteady  transonic  flows.  However,  at  the  same  time  there 
is  a  shortage  of  experimental  results  [1].  The  capabilities  of  ground 
test  facilities  have  actually  been  surpassed  by  analytical  and  computa¬ 
tional  methods  in  many  areas  of  unsteady  aerodynamics.  One  reason  is  the 
inherent  problem  of  obtaining  detailed  measurements  about  a  moving  model. 
However,  a  large  class  of  unsteady  flow  experiments  can  be  conducted  with 
stationary  models  when  a  controlled  oscillatory  flow  is  generated  in  the 
test  section.  These  so  called  gust  tunnels  have  demonstrated  the 
potential  of  providing  valuable  aerodynamic  data  [2]. 

Unfortunately,  with  few  exceptions  [3],  all  existing  gust  tunnels  are 
small,  low  speed  facilities.  Recently,  however,  a  new  concept  termed  the 
"Ball  Wall"  has  been  used  to  introduce  oscillations  in  the  transverse 
velocity  over  a  wide  range  of  speeds  [4].  The  new  concept  utilizes 
a  temporal  variation  of  the  test  section  wall  porosity  to  produce  a  fre¬ 
quency  range  of  flow  perturbations  within  the  test  section  that  is  unique 
relative  to  other  gust  tunnels.  Measureable  flow  angularity  oscillations 
were  generated  at  frequencies  from  quasi-steady  to  50  Hz  in  tunnel  {4T)  at 
the  Arnold  Engineering  Development  Center  [4]. 

Having  demonstrated  the  experimental  capability  of  producing  a  con¬ 
trolled  oscillatory  flow  from  subsonic  incompressible  to  the  transonic 
flow  regime,  a  serious  question  arose  concerning  the  dynamic  similarity 
due  to  an  oscillating  flow  as  opposed  to  an  oscillating  model  of  equal 
frequency  and  amplitude.  Specifically,  the  question  is:  is  the  velocity 
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distribution  In  the  neighborhood  of  the  model  and  the  resulting  normal  and 
shear  stress  distribution  equivalent  when  the  free  stream  contains  a 
periodic  transverse  oscillation  which  is  equal  in  magnitude  to  transverse 
model  oscillation? 

It  is  well  known  [5]  that  for  compressible  flow  the  two  problems  of 
fixed  bodies  in  oscillating  flows  and  oscillation  bodies  in  steady  uniform 
flows  are  not  equivalent.  The  discrepancy  is  due  to  the  inertia  term  in¬ 
troduced  in  the  boundary  layer  equations  when  transforming  from  the  fixed 
to  the  oscillating  coordinate  system.  Thus,  the  boundary  layer  equations 
are  different, and  as  expected  the  near  field  velocity  distribution  is  not 
equivalent  for  the  two  cases. 

A  common  misconception  is  that  incompressible  flows  the  two  prob¬ 
lems  are  always  equivalent,  although  for  certain  in-plane  oscillations 
equivalent  surface  stress  distributions  are  produced.  One  example  of 
equivalent  transverse  oscillations  is  the  two-dimensional  stagnation  point 
flow  generated  by  an  infinite  flat  plate.  A  uniform  flow,  U^,  approaching 
a  plane  wall  which  oscillates  with  velocity,  Vp^,  given  by  Vp  =  cos  ut 
produces  the  equivalent  wall  shear  stress  distribution  as  that  of  a  fixed 
wall  which  sees  a  uniform  flow  with  an  equivalent  transverse  component 
which  is  It  out  of  phase.  These  two  problems  were  first  investigated  by 
Glauert  [6]  and  the  kinematic  and  geometric  similarity  which  produces  the 
equivalent  dynamic  effects  is  physically  obvious.  In  fact  a  simple 
coordinate  transformation  from  one  fixed  in  space  to  one  fixed  to  the 
stagnation  point  demonstrates  the  equivalence  [6]. 

Now  consider  the  two  distinct  problems  of  a  right  circular  cylinder 
undergoing  small  transverse  oscillations  in  a  steady  uniform  flow  and  what 
appears  to  be  the  similar  problem  of  a  stationary  cylinder  subjected  to  a 


steady  uniform  flow  which  contains  the  equivalent  transverse  oscillations, 
i.e.  equal  frequency  and  magnitude.  For  the  idealized  situation  all 
points  in  the  flow  field  have  the  same  vertical  oscillation,  which  is 
different  from  the  so-called  gust  problem  where  a  periodic  component  is 
convected  down  stream.  The  gust  condition  may  be  a  more  realistic  model 
for  the  conditions  created  in  some  wind  tunnels,  however  it  is  obvious 
that  the  flow  field  for  this  condition  is  not  even  kinematically  similar 
to  the  case  of  an  oscillating  cylinder. 

If  the  potential  flow  field  due  to  a  steady  uniform  form,  U  ,  plus  a 

00 

small  transverse  oscillating  component,  given  by  cos  J-  is  solved,  the 
resulting  surface  velocity  for  an(r,e)  coordinate  system  fixed  to  the 
cylinder  is 


Vg(R,e,t)  =  -2U„  Sin  e  -  Coswt  Cos  e. 

The  stagnation  point  is  then  displaced  through  maximum  angle  given  by 

^S.P.  “  JT  ' 
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Now  consider  the  corresponding  potential  flow  problem  of  a  uniform 
steady  flow  approaching  a  cylinder  undergoing  a  vertical  oscillation  of 
identical  frequency  and  amplitude  but  with  a  180°  phase  difference.  Solv¬ 
ing  this  problem  by  an  analysis  similar  to  that  of  G.  I.  Taylor  [7],  the 
resulting  surface  velocity  is  given  by 


Vg(R,0,t)  =  2U^  Sin  e  -  Cosut  Cos  e. 

It  is  seen  that  the  unsteady  part  of  v  (R,e,t)  is  different  and  in  fact 
the  maximum  stagnation  point  excursion  is  exactly  one-half  that  of  the 
first  case.  Jt  then  follows  that  the  surface  shear  stress  distribution 
will  also  be  different,  and  in  the  linearized  case  of  small  transverse 
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oscillations  is  one-half  that  of  the  first  case.  The  reason  for  the 
factor  of  two  discrepancy  is  due  to  the  physical  difference  in  the  two 
flow  fields.  In  the  first  problem,  the  entire  flow  field  oscillates  in 
time  and  the  presence  of  the  cylinder  is  essentially  ignored.  There  is  no 
circulation  and  the  flow  is  symmetric  about  the  dividing  streamline,  a 
condition  which  closely  approximates  the  quasi-steady  viscous  problem.  In 
the  second  case  the  flow  field  is  basically  steady,  and  the  unsteady 
effects  are  limited  to  a  small  region  around  the  cylinder,  known  as  the 
added  mass  region.  An  observer  fixed  on  the  cylinder  sees  a  negative 
velocity,  V^,  at  infinity,  however  as  the  particles  approach  the  cylinder, 
they  experience  a  slowly  increasing  positive  velocity  field.  The  net 
effect  is  to  produce  an  asymmetric  flow  and  a  stagnation  point  deflection 
of  one-half  that  of  the  first  case. 

From  this  example,  even  when  all  points  in  the  flow  field  contain  a 
spatially  uniform  and  instantaneously  oscillating  transverse  component,  it 
is  seen  that  there  is  a  difference  in  the  near  field  fluid  mechanics.  The 


differences  can  become  more  pronounced, when  in  fact  the  spatial  distribu¬ 
tion  is  not  quite  uniform.  This  is  in  fact  what  occurs  near  the  beginning 
of  the  porous  wall  region,  especially  when  the  freestream  velocity  is 
large  [4]  and  the  length  of  the  porous  wall  is  of  the  same  order  as  the 


model . 

Having  gained  some  physical  insight  into  the  subtle  differences  of 
the  two  unsteady  problems,  i.e.  oscillation  flow  fields  vs.  oscillating 
bodies,  during  the  AFOSR  Summer  Research  it  was  decided  to 

study  the  unsteady  fluid  dynamics  around  a  restricted  class  of  airfoils. 
Specifically,  the  investigation  is  limited  to  airfoils  having  small  thick- 
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the  subsequent  numerical  integration  of  the  boundary  layer  equations. 
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2.0  MATHEMICAL  FORMULATION 


The  basic  problem  considered  in  this  investigation  is  an  unsteady 
laminar  incompressible  flow  over  a  two-dimensional  thin  airfoil.  The 
concept  of  a  thin  airfoil  as  used  in  this  analysis  implies  that  t/b  <<  1 
as  shown  in  Figure  1.  This  implies  that  the  mean  chord  of  the  airfoil 
lies  along  the  X  axis  to  within  0[  e],  where  e  =  a/b  and  a/t  =  0[1].  This 
implies  that  disturbances  to  the  free  stream  flow  generated  by  the  air¬ 
foil  are  0[e].  This  formulation  has  distinct  mathematical  advantages 
which  will  be  explored  in  the  following  analysis.  First,  it  allows  for 
an  analytical  solution  to  the  potential  flow  field  to  within  0[e  ].  And 
secondly, it  permits  a  very  accurate  and  stable  hybrid  numerical  scheme  to 
be  employed  for  the  solution  of  the  boundary  layer  equations.  Although 
the  thin  airfoil  approximation  limits  the  physical  application  of  the 
analysis,  its  simplicity  allows  for  both  physical  insight  into  the  un¬ 
steady  efforts,  and  provides  a  valuable  base  line  analysis  for  comparison 
to  subsequent  studies  which  can  include  airfoil  thickness,  camber,  and 
angle  of  attack. 

Specifically,  the  purpose  of  this  analysis  is  to  compare  the  effects 
of  two  distinct  unsteady  flows  under  the  thin  airfoil  approximation. 

First,  a  uniform  free  stream  approaches  the  airfoil  which  is  undergoing  a 
periodic  transverse  oscillation  with  respect  to  the  x-y  axis  of  amplitude 
a  and  frequency  a,  as  shown  in  Figure  2a.  And  secondly,  a  uniform  flow 
plus  a  periodic  transverse  component  approach  the  same  airfoil  which  is 
now  fixed  with  respect  to  the  X-Y  axis.  This  part  of  the  analysis  is  suf¬ 
ficiently  general  so  that  any  spatial  distribution  in  the  transverse  com¬ 
ponent  capable  of  beign  represented  by  a  Fourier  series  can  be  consider- 
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ed.  For  both  conditions  the  amplitude  of  the  oscillation  is  "small", 
i.e.,  a/b  «  1.  This  implies  that  both  the  airfoil  geometry  and  the 
unsteady  effect  will  appear  first  in  the  0[e]  equations. 

2.1  Potential  Flow  Formulation 

In  addition  to  the  previous  restrictions,  the  fluid  will  be  assumed 
inviscid  in  this  section.  This  is  reasonable  under  thin  airfoil  theory, 
since  without  separation  the  viscous  effects  are  confined  to  the  boundary 
layer  and  the  effect  of  viscosity  can  be  incorporated  by  a  phenomenologi¬ 
cal  rule  that  the  velocity  at  the  trailing  edge  must  remain  finite  and  tangent 
to  the  airfoil  surface.  Finally,  in  addition  to  the  inviscid  assumption, 
the  fluid  will  be  taken  as  piezotropic,  i.e.,  o  =  p(p)  only. 

The  equations  will  be  formulated  in  terms  of  the  acceleration 
potential,  <>,  of  the  flow  field.  This  approach  is  selected  because  if  is 
continuous  throughout  the  idealized  wake  in  thin  airfoil  theory  whereas 
the  usual  velocity  potential  is  discontinuous.  A  detailed  discussion 
of  the  derivation  of  the  equations  can  be  found  in  Ref.  [8].  In  this  sec¬ 
tion,  we  simply  note  that  the  relavent  equations  are  Euler’s  equation 
written  as 


Ot 


and,  the  acceleration  potential  definition  is  given  by 


DV 

"  m 


(1) 


(2) 


Substituting  equation  (2)  into  (1)  and  integrating  with  respect  to  space 
yields 
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■  ^ 
p 


constant 


Further  simplifications  plus  a  discussion  of  the  boundary  conditions  will 
be  made  in  Section  2.3. 

2.2.  Boundary  Layer  Formulation 

The  two  dimensional*  incompressible,  unsteady,  laminar  boundary  layer 
equations  in  terms  of  the  stream  function  are 


2  2  2  3 

3  ij;  ^  3  Ip  3  iji  _  3'P  _ i 

3t3y  Ty  3x3y  lx  ~  p  3x  ''  ^^3 


where  ^  —  and  v  h  -  . 

The  unsteady  pressure  gradient  is  given  by 


,  3U  3U 

I  i£  =  _£  +  u  — 

p  3X  3t  e  3X 


Adopting  the  usual  body  coordinate  system  shown  in  Figure  3,  the 
boundary  conditions  become 


^  ^  =  0  at  y  =  0 

3y  3X 


-  Ug(x,t)  as  y 


(6a, b) 


In  equation  (6c),  U^(x,t)  represents  the  velocity  at  the  edge  of  the 
boundary  layer  which  will  be  determined  from  the  potential  flow  solution. 
2.3.  Nondimensionalization  and  Perturbation  Equations 

Outside  the  leading  edge  region,  two  important  length  scales 
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characterize  the  unsteady  boundary  layer.  They  are  6^  and  Sg,  where 
6-  -  (v/uj)  '  and  6-  ~  (vx/U  )  .  If  no  steady  flow  is  present, 

i  D  “ 

the  fluctuating  part  of  the  flow  extends  a  distaance 6^  away  from  the  sur¬ 
face.  When  a  steady  flow  is  superimposed,  the  two  layers  interact  and 
the  resulting  flow  can  be  quite  complicated  when  the  two  layers  are  of 
comparable  magnitude,  i.e.  —  =  0[1].  However,  under  certain  restric- 

Utf} 

tions  the  mathematical  analysis  can  be  dramatically  simplified.  This 
occurs  if  we  restrict  the  magnitude  of  the  transverse  oscillating  com¬ 
ponent,  Vq,  such  that  V^/U^  <<  1.  This  does  not  limit  the  analysis 
since,  by  definition  of  the  problem,  is  less  than  about  10%  of  the 
free  stream  velocity  due  to  physical  limitations  of  the  wind  tunnel. 
Referring  to  Figure  2,  we  see  that  -  i*)a,  thus  uia/U^  <<  1.  This  can  be 
written  as  (dib/U  )e  «  1,  where  e  h  a/b.  This  implies  i*)b/U  =  0[1]  or 
equivalently  wx/U^  =  0[1]  except  near  the  leading  edge  of  the  airfoil. 
However  by  definition,  the  boundary  layer  analysis  excludes  this  region, 
thus  we  conclude  ojx/U  =  0[1]. 

OD 

The  parameter  has  a  special  significant,  consequently  we  let 

00 

a  =  2bw/U^,  and  a  =  0  [1]  implies  that  the  two  layers  are  of  comparable 
magnitude,  and  either  6^  or  6g  are  acceptable  length  scales.  However, 
since  e  ->•  0  implies  the  oscillatory  component  vanishes,  Sg  is  more 
appropriate.  Accordingly  the  following  nondimensional  variables  and 
parameters  are  selected  for  the  boundary  layer  phase  of  the  analysis. 


flo 
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R 


e 


i>  5 


(2bvU  ) 


T7? 


2bu>  3 

^  -  '5  ' 

do 


(7) 


Substituting  equations  (7)  into  equation  (4)  yields 


2-  -  2-  -  2-  3- 

3  i(/  ^  3  -  _3iti3<;_3ip__  3p 

TrTn  3n  StSn  "sT  .  2  T~5  ac 


(8) 


where 


-  ii  = 

3c 


a 


+  u 


e  3c 


(9) 


with  the  following  boundary  conditions. 


it 

3r, 

=  it  =  0 

3C 

at 

n  =  0 

(10a, b) 

it 

3n 

^  ^e 

as 

T)  ao 

(10c) 

From  the  physical  description  of  the  problem,  it  is  reasonable  to 
expect  that  both  the  edge  velocity,  U^,  and  the  pressure  gradient  calcu¬ 
lated  from  potential  flow  theory  can  be  expressed  as 

Og(c,T)  =  Og  +  eU,(c,T)  +  0[e^] 

and  _  « 

P(c.t)  =  Pq  +  eP^(c,T)  +  0[c^] 

It  then  follows  that  the  stream  function  can  be  written  as 

-  -  -  2 

<<(c,n,T;a,6)  =  4'Q(c»n;a)  +  ei|(^  (c,n,-t;a)  +  0[e  ]  (13) 


(11) 

(12) 


11 


Substituting  equation  (13)  into  the 

stream  function  equation  yields  the 

following  two  equations. 

(14) 

TiT’  3C3n  3c  ,  2 

0  n 

'IT'  ■ 

and 

2  2-  * 

3  i|>i  ^  '1)1  _ 

a - -  +  i  - 

3T3n  o,n  3;an  o,5n  3n 

34;, 

_  4,  0  ^ 

0,c  .  2  0,nn  3c 

3n 

3^4^  SPj 

3ri^ 

(15) 

In  equation  (15)  the  notation 

0,  n 

E  3  4>Q/3n  is  adopted  for 

simplicity. 

The  corresponding  set  of  boundary  conditions  become 

_2  =  —2  =  0 

3ri  3n 

at  n  =  0 

(16a, b) 

— 2  -►  u 

3n  0 

as  n  ® 

(16c) 

and 

3'I'l  3^^ 

ItT  '3c  '  ° 

ai*. 

at  n  =  0 

{17a, b) 

as  n  ♦ 

(17c) 

In  the  potential  flow  region,  the  characteristic  length  scale  for 
both  X  and  Y  is  simply  one-half  the  airfoil  chord,  b.  The  velocity  scale 
is  the  free  stream  velocity,  U.  The  time  and  acceleration  potential  are 
nondimensionalized  by  b/  and  respectively.  All  other  nondimensional 
variables  are  unchanged.  The  linearized  potential  flow  equations  can  now 
be  easily  formulated  by  again  expanding  the  nondimensional  pressure,  den¬ 
sity  and  acceleration  potential  in  terms  of  t.  Substituting  these 
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expansions  into  equation  (8)  yields 

^-1  dP  dP 


$  +  e 

0 


O  #  ’  ur 

5  +  0[c2]  =  =  [-° 

•  J 


"  Cl 


-  c  J-  +  oCt'^]] 


(18) 


Thus, 


=  constant,  for  P^  =  constant;  and  41^  =  -  P-j/pQ  • 

The  linearized  field  equations  can  also  be  easily  obtained.  First, 
Eulers  equation  becomes 


au 


at 


I-  ^ 

^0  aX 


av 


and 


at 


pq  aY 


(19a, b) 


where  the  convective  acceleration  terms  are  OCe'^].  Similarly,  the 
continuity  equation  becomes 


;  (ii  .  iv,  ,  !!i  .  0 

°  aX  aY  at 


(20) 


Equations  (19)  and  (20)  can  be  combined  to  yield 


a^P 

,2p 

'  a"  at"^ 


(21) 


For  an  incompressible  fluid  the  speed  of  sound,  a,  approaches 
infinity  and  the  linearized  equation  becomes 


=  0 


(22) 


where  the  substitution  5^  =  P^/Pq  was  made  in  equation  (22). 
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Finally,  the  acceleration  components  can  be  linearized  as  follows: 
From  equation  (2)  for  two  dimensional  flow  we  obtain 

ii  =  and  ii  = 

3X  3t  sY  Dt 

Substituting  the  asymptotic  expansion  for  $  and  noting  that  u 
£U^(X,Y,t)  +  0(e  ),  and  equating  equal  powers  of  f  yields 

3$^  _  3Ui  3U-J 

3X  °  3X  3t 

Similarly,  for  3(t)^/3Y,  we  obtain 

30^  _  3V^  3V^ 

sY  °  3X  3t 

Having  formulated  the  general  linearized  unsteady  equations  the  solution 
for  the  potential  flow  region  and  the  boundary  layer  will  be  discussed 
separately. 


3.0  POTENTIAL  FLOW  SOLUTION 


Within  the  framework  of  linearized  theory,  individual  solutions  may 
be  superimposed  to  generate  the  complete  solution.  Thus, the  flow  field  of 
an  oscillating  airfoil  with  "small"  but  finite  thickness  and  camber  can  be 
obtained  by  the  superposition  of  an  unsteady  airfoil  of  zero  thickness  and 
camber  and  a  steady-state  solution  for  an  airfoil  of  small  thickness  and 
camber.  However,  in  this  analysis  the  characteristics  of  the  two  distinct 
unsteady  effects,  i.e.,  oscillating  airfoil  vs  oscillating  free  stream, 
were  of  primary  interest.  For  this  reason,  as  a  first  approximation,  it 
was  decided  to  look  at  the  steady  state  component  for  the  airfoil  with 
zero  thickness  and  camber.  This  approximation  permits  closed  form 
analytical  solutions  to  the  potential  flow.  This  provides  valuable 
physical  insight  into  the  role  of  the  bound  vorticity  and  the  added  mass 
effect.  In  addition,  it  provides  a  reference  point  for  further  investiga¬ 
tions  which  can  include  both  small  but  finite  thickness  and  camber 
effects. 

3.1  Thin  Airfoil  Solution 


Under  these  restrictions,  the  equation  describing  the  unsteady  ef¬ 
fects  can  be  written  below,  where  here  and  in  subsequent  equations  the  over¬ 
bar  for  the  nondimensional  terms  is  dropped. 

2  2 

!!i  +  !!i  -  0  <26) 

3X2  3y2“ 

The  steady  state  velocity  is  independent  of  X,  and  the  linearized  accel¬ 
eration  components  are  given  by  equations  (24)  and  (25). 

Refering  to  Figure  1,  we  can  write  the  boundary  conditions  at  the 
airfoil  surface  as 
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)  _  Oh  _  3h  ,,  ah 

'''rs  ■  TJt  '  m  ^  %  37 


(27  ) 


and 


fri 


D^h  _  3^h  y..  3^h  ,,2  3^h 


(28) 


where  h(X,t)  describes  the  position  of  the  airfoil.  In  addition  to  these 
conditions  we  require  and  -  0  as  (X  +  Y  )  plus  41^  must  be 

finite  at  the  trailing  edge.  This  last  condition  is  known  as  the  Kutta- 
Joukowski  condition.  Unfortunately,  due  to  the  conformal  transformation 
the  acceleration  potential  will  not  be  finite  at  the  leading  edge.  This 
singularity  and  its  subsequent  elimination  will  be  discussed  in  Section 
3.2. 

It  is  now  that  the  powerful  mathematical  techniques  of  potential 
theory  can  be  exploited.  It  is  well  known  that  a  two  dimensional  airfoil 
can  be  treated  by  conformal  mapping.  Since  the  linearized  boundary  condi¬ 
tions  are  expanded  about <  X  <  1  gpd  Y  =  0,  we  can  map  this  into  a  circle 
by 


Z  =  '  (Z*  *  4) 

2  z 

Thus  the  boundary  conditions  in  the  physical  Z  plane  can  be  written  in  the 


conformal  Z  plane  (see  figure  4)  as 


'1  -1  ^'*’1 

(jjp)s  (r=l,  e=Cos  X,t)  =  (-^)5  (X,o,t)  Sine 


(29) 


In  addition,  the  condition  at  infinity  can  be  written  as 


lim  Real  {wl  =  constant 

1 2* I  *  ” 


(30) 
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where  w  is  the  complex  acceleration  potential  given  by 


w  =  (ti^  +  i’l'^  =  f(X  4  iY) 

* 

In  the  Z  plane,  w  has  the  form 


(31) 


w(Z*)  =  i  -5P^  +  i  ^  (32) 

Z  +1  Z 

where  A  and  B  are  determined  from  the  prescribed  motion  of  the  airfoil  by 
applying  the  boundary  and  Kutta  conditions. 

3.1.1  Vertical  Oscillation  of  Airfoil 

For  an  airfoil  undergoing  vertical  translation  oscillations  described 
by  h  =  a  Cos  ut,  the  solution  to  the  acceleration  potential  is  given  by 

—  =  {-ik  C(k)  Tan(^°^-  -)  +  k^  Sin(Cos‘^X) }e  ^  (33) 

Pq  2 

where  X  =  2c  -1  (see  Figure  3),  K  =  -^  =  ^  and  C(k)  =  F(k)  +  iG(k).  The 
function  C(k)  is  known  as  Theodorsen's  function  and  its  value  along  with 
the  details  of  the  solution  can  be  found  in  Ref.  [9].  From  the  defini¬ 
tion  of  we  see  that  the  unsteady  pressure  gradient,  which  is  of 
interest  in  the  boundary  layer  analysis  is  simply 


(34) 


Performing  the  indicated  differentiation,  and  after  lenghty  simpli¬ 
fications  we  obtain  the  unsteady  pressure  gradient  in  boundary  layer  co¬ 
ordinates  given  by 


1 


KaG  -  -  1t.F)e” 


The  corresponding  "edge"  velocity  can  then  be  found  from  the  unsteady 


Bernoulli  equation  in  nondimensional  boundary  layer  coordinates  given  by 
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Equation  (9). 

At  this  point,  the  potential  analysis  would  be  complete  for  the 

oscillating  airfoil  if  it  were  not  for  the  leading  edge,  i.e.  t  ♦  0, 

singularity  which  is  apparent  from  equation  (35), 

3.1.2  Vertical  Oscillation  of  Free  Stream 

The  second  unsteady  flow  field  in  this  investigation  is  that  of  a 

uniform  flow  in  the  X  direction  plus  a  small  periodic  component  in  the  Y 

direction  as  shown  in  Figure  2b.  With  coordinates  fixed  on  the  airfoil, 

the  effect  of  this  oscillation  in  dimensional  form  can  be  expressed  by 
i 

Y  =  h(X)e  .  The  boundary  condition  on  the  airfoil  surface  given  by  equa¬ 
tion  (26)  becomes 

.a . 

1-xt 

In  nondimensional  form  we  have  Y  =  h(X)  e  ,  thus 


<''1>S 


(i  ^  h  +  ^)e 


(37) 


The  unsteady  transverse  velocity  may  contain  a  spatial  nonuniformity 
due  to  the  physical  limitations  of  the  wind  tunnel.  In  fact  it  is  difficult 
to  create  a  uniform  transverse  oscillating  flow  at  the  wind  tunnel  centerline 
Typical  transverse  spatial  distributions  are  shown  in  Figures  6  and  7  and  the 
physical  characteristics  are  discussed  in  Section  4.3.  For  computational 
purposes,  h(X)  is  taken  to  be  CosX,  which  closely  approximates  the  distribu¬ 
tion  shown  in  Figure  6  and  represents  one  of  the  more  nonuniform  conditions. 
For  this  distribution,  the  "upwash"  on  the  airfoil  becomes 


(Vi)s  =  (i|  CosX  -  SinX)e 


i^t 


(38) 


Again,  omitting  the  details  which  can  be  found  in  Ref.  [8]  or  [9],  we  obtain 


T 


the  unsteady  pressure  distribution  given  by 


,Cos"’x 


-1. 


P,  =  (2a^  Tan(^”«  --)  +  4  I  a„  Sin(nCos  'X))e 
'0  2  n 


i|t 


(39) 


where 


%  ‘  <=<''"^'’0  *  I’ll  -  fi 

»n  '  S  ".-1  *  'n  -  n  '■-,*1  •  "  i  ’ 

The  coefficients  are  defined  in  terms  of  the  "gust"  velocity. 

the  relation  (see  Figure  4)  e  =  cos  the  coefficients  are 

n 

P.  =  -  f  V, (e ,t)  do 

0  n  J  I 


P 


n 


0 

n 

-  ~  [  V,  (e  ,t)  Cos{n  e)  de  , 

n  j  I 


n  >  1 


(40) 

(41) 

Using 


(42) 

(43) 


0 

From  this,  the  unsteady  pressure  gradient  in  boundary  layer  coordinates 
(i;,ri,T)  becomes 

JT  "  ^  -37?-- - 177  ^  a  [n  Cos(n  Cos'^  (2c-l)])e^^ 

Equation  (44)  is  the  counterpart  to  equation  (39)  for  an  oscillating 
airfoil  in  a  steady  uniform  flow.  In  fact  equation  (44)  does  reduce  to 
equation  (39)  if  h(x)  a  which  implies  (v^l^  +  i^  Then  P^  =  -i^ 

and  Pp  =  0  for  n  >  1  which  implies  f  ®1  " 

Since  the  pressure  gradient  given  by  equation  (44)  can  represent 
either  of  the  two  distinct  unsteady  flows,  it  was  used  in  the  boundary  con¬ 
dition  part  of  the  boundary  layer  code  to  calculate  the  edge  velocity 
Ug(c,T).  A  ten  point  Gausian  quadrative  scheme  was  used,  and  the  first 
twenty  coefficients  were  calculated.  The  truncation  accuracy  of  the 

series  including  twenty  coefficients  insured  overall  relative  accuracy  of 
*8 

10  which  is  more  than  sufficient  for  the  analysis. 
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As  expected,  the  unsteady  pressure  gradient  due  to  a  periodic  trans¬ 
verse  oscillation  convected  over  a  stationary  airfoil  contains  a  leading 
edge  singularity.  This  mathematical  singularity  must  be  avoided  or  re¬ 
moved  in  some  logical  manner  so  that  the  boundary  layer  equations  can  be 
integrated.  The  method  selected  to  accomplish  this  is  discussed  in  the 
next  section. 

3.2  Leading  Edge  Solution 

_2 

Due  to  the  conformal  mapping,  the  pressure  gradient  has  a  ^  singu- 

-1/2 

larity  at  the  leading  edge  and  a  ^  singularity  at  the  trailing  edge. 

_2 

The  trailing  edge  mathematical  singularity  is  convergent,  however,  the  ^ 
is  not.  One  obvious  way  of  avoiding  this  difficulty,  which  is  suggested  in 
Ref.  [7]  is  to  simply  start  the  integration  at  a  small  distance  from  the 
leading  edge  at,  say  4%  chord.  Although  this  does  not  significantly  effect 
the  unsteady  lift  calculations,  it  does  produce  significant  changes  in  the 
boundary  layer  calculations  and  hence  the  unsteady  drag.  For  this  reason, 
the  following  method  was  selected  for  starting  the  boundary  layer  calcula¬ 
tions:  It  is  well  known  that  the  leading  edge  region  of  a  two-dimensional 
airfoil  can  be  approximated  by  the  stagnation  flow  around  a  right  circular 
cylinder.  Consequently,  the  flow  field  at  the  stagnation  point  of  a 

cylinder  oscilliating  with  the  same  amplitude  and  frequency  as  the  airfoil 
can  be  used  as  initial  values  for  both  and  dP/dx  in  the  boundary  layer 
calcuations. 

For  a  cylinder  oscillating  vertically  in  an  inviscid  uniform  flow,  the 
appropriate  field  equation  in  terms  of  the  stream  function  is  simply 

2 

=  0  (45) 

The  appropriate  boundary  condition  on  the  surface  in  (r,o)  coordinates  is 


~  fft  '  Sin  e 
r  36  ■'o 


(46) 


where  represents  its  velocity  with  respect  to  (X,Y)  as  shown  in  Figure  5. 

Following  the  solution  first  described  by  G.  I.  Taylor  [7],  the  cir¬ 
cumferential  velocity  at  r  =  R  is  given  by 


V  (R,e)  =  -y  Cos  e  -  2U  Sin  e 
6  0 


(47) 


Letting  yQ  =  a  cos  ait,  we  can  easily  write  the  edge  velocity,  U^,  in  bound¬ 


ary  layer  coordinates  as 


Ug(x,t)  =  2U^  Sin  ^  +  3“  Cos  ^  Sin  a)t 


(48) 


At  the  stagnation  point  sin  ^  and  equation  (48)  becomes 

'r 


U„(c.t)  =  4mc  -  i-^ee^^ 


149) 


where  Equation  (49)  is  nondimensional  and  m  =  b/R. 

The  unsteady  pressure  gradient  can  now  be  found  from  equation  (9)  by 


expanding  both  P  and  in  terms  of  e ,  which  yields 
3P,  3U,  3U,  3U„ 

-  — L  =  +  u  -I  +  u,  -2 

3i;  3t  0  3t  1  3c 


(50) 


Noting  from  equation  (49)  that  =  4m^  and  =  i  “  and  substituting 
into  equation  (50)  yields 


TT  =  (f-  • 


(51) 


The  potential  flow  field  is  now  completely  specified  by  the  pressure 
gradients  given  by  equations  (35)  and  (44)  plus  the  stagnation  point 
pressure  gradient  given  by  equation  (51).  It  is  interesting  to  compare 
equation  (35)  in  the  limit  as  i;  0.  Looking  at  the  real  part  of  sP^/ac 
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which  represents  the  added  mass  effect,  i.e.  neglecting  the  bound  vorticity 
component,  we  have 

dp,  2 

lim  Real  ~ 

c^O 

This  is  precisely  the  real  part  for  the  cylinder  given  by  equation  (51). 
This  is  to  be  expected,  since  the  cylinder  has  no  circulation  and  we  see 
that  the  singularity  in  equation  (35)  is  due  to  the  bound  vorticity 
component. 
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4.0  BOUNDARY  LAYER  RESULTS 


4.1  Thin  Airfoil  Solution 

The  solution  to  the  boundary  layer  outside  the  leading  edge  region 
through  0[e]  is  given  by  integrating  equations  (14)  and  (15)  with  the 
appropriate  boundary  conditions  and  using  the  results  of  Section  3.1. 
Equation  (14)  represents  the  steady  flow  contribution  for  an  airfoil  of  zero 
thickness  and  hence  zero  pressure  gradient.  This  equation  can  be  rewritten 
in  the  usual  form  by  introducing  the  Blasius  similarity  variables  given  by 


(vx/U  ) 


m 


and 


’^B  " 


'i'(x,y) 


vxU  ) 


m 


(vxU  ) 


TTT 


From  these  two  relations,  it  follows  that 


n 


and  i<^(c,n)  = 


(53a, b) 


Substituting  equations  (53a, b)  into  Equation  (14)  yields 


1  '  ‘ 
—  f  f 
2  B^B 


=  0 


(54. ) 


with  boundary  conditions  given  by 
fg(o)  =  fg(o)  =  0 


(55a, b) 


t 

fglng-)  >  1 

The  result,  as  expected,  is  the  well  known  Blasius  solution. 


(55c) 
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Turning  now  to  the  unsteady  component  and  using  u^ 
Vq  ^  equation  (15)  can  be  written  as 

Zii, 


'oVq^  3r  and 


3^^  3^4-1  “''1 

_ 1  ,  ,,  _ L  -  M  -  +  V 

3r,3T  0  3;3n  o,n  3n  0 


(66) 


an 


,3 


-  u 


a  ij/l  a  ■ 


o,n  a; 


3r,' 


ap^ 

TT 


where 


u_  =  fr 


and 


1 

T7? 


-  f. 


(57a, b) 


From  equations  (57)  we  see  that  in  addition  to  the  singularity  in  the  pres¬ 
sure  gradient  at  ?  =  0,  the  coefficients  also  contain  singularities.  This 
presents  no  additional  difficulty,  since  at  =  0  the  corresponding  stagna¬ 
tion  point  boundary  layer  equations  along  with  the  pressure  gradient  given 
by  equation  (51)  are  used  to  start  the  integration. 

The  temporal  variation  in  equation  (56)  can  be  separated  by  the  complex 
exponential  time  factor.  This  is  due  to  the  linearity  of  the  equation  which 
is  a  consequence  of  the  small  amplitude  transverse  oscillation.  Represent¬ 
ing  the  stream  function  and  the  pressure  gradient  in  the  complex  plane  as 

';<^(c,ri,T)  =  (c.nle^"' 


and 

P^(;,t)  =  P^(c)e^" 

and  substituting  these  into  equation  (56)  yields 


la 


ay,  a^y, 

— L  +  jj  - L  -  V 

an  0  acan  o,n 


ay^  a^y^ 

^  aT" 


ay 


o,n  ac 


3 

3'^y, 


an" 


IT 


(59) 


(50) 


This  equation  could  be  integrated  by  a  straightforward  finite  dif¬ 
ference  method,  however  a  more  flexible  hybrid  scheme  is  adopted.  First, 
the  stream  function  and  the  pressure  gradient  are  decomposed  into  their  real 
and  inaginary  parts,  i.e.,  4'  =  +  iy_  Then,  only  the  ■'.-direction  is 

finite  differenced  according  to 


+  0(AC) 


(61) 


Substituting  these  results  into  equation  (60)  yields  the  following  two  equa- 
rions  in  n: 


0  r 


0,  n 


^)y'  + 

Ac  r 


0,  n 
Ac 


+  ay^ 


dp„  n 


r  Q  ui  .  0 » 

AC  r  AC 


n-1 


y  } 

r 


(62) 


''o^ 


+  (v  -  — )y-  + 

0,n  AC  1 


0,n 

Ac 


-  ay 


.  dp  n 
r  ■ 


AC  I 


U 

o,n 

Ac 


(63) 


with  the  following  boundary  conditions: 


o 

II 

at 

n  ~  0 

(64a, b) 

U^{c) 

as 

n  ->  m 

(64c) 

y^.  =  0 

at 

n  =  0 

(65a, b) 

U.(c) 

as 

T)  ^  ® 

(65c) 
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Equations  (62)  through  (65)  constitutes  a  well  posed  coupled  boundary 
value  problem.  Its  solution  is  obtained  by  using  a  general  purpose  spline 
collocation  code.  This  particular  integration  scheme  was  selected  because  it 
is  more  versatile,  accurate  and  stable  than  a  full  finite  difference 
approach  when  solving  parabolic  equations.  The  accuracy  is  controlled  by 
the  step  size,  since  the  spline  collocation  integration  in  the  n  direc¬ 
tion  is  essentially  exact  (relative  error  tolerance  on  the  function  can  be 
set  at  10  ^  without  convergence  difficulties).  In  addition,  as  A;  -►  0  the 
finite  difference  characteristics  approach  the  continuum  characteristics 
which  are  normal  to  the  caxis.  A  final  justification  for  adopting  this 
scheme  over  full  finite  difference  methods  lies  in  its  ability  to  easily 
incorporate  moving  boundaries.  Recall  that  in  the  present  formulation,  the 
boundary  conditions  are  applied  to  the  mean  position,  i.e.  Y  =  0,  of  the 
oscillating  airfoil.  However,  it  would  be  of  interest  in  future  investiga¬ 
tions  to  study  the  effect  of  applying  the  boundary  condition  to  the  moving 
surface  while  holding  all  other  parameters  constant.  This  would  then 
require  only  slight  modifications  to  the  boundary  layer  code  whereas  for 
full  finite  difference  methods  moving  boundaries  would  introduce  major 
complications. 

Before  discussing  the  results  of  the  numerical  solution  of  equations 
(62)  through  (65)  the  equations  valid  in  the  neighborhood  of  c =  0  wi 1 1  be 
developed.  The  solution  of  these  stagnation  point  equations  will  provide 
the  starting  profile  for  the  boundary  layer  calculations. 

4.2  Leading  Edge  Solution. 

The  flow  field  near  the  leading  edge  is  described  by  equations  (14) 
through  (17)  along  with  the  potential  flow  results  of  Section  3.2.  The  same 
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Now  writing  and  as  and  P^(de^"  we  obtain 

the  equation  for  the  stagnation  point  region  given  by 


ia  ^  +  4n,f  •  ^  .  2m 


Notice  that  equation  (71)  is  an  ordinary  differential  equation  whereas  equa¬ 
tion  (60)  is  a  partial  differential  equation.  This  is  due  to  the  fact  that 
=  i  2  which  is  independent  of  c.  Equation  (71)  now  represents  a 
boundary  value  problem  where  the  boundary  conditions  are 


$l(0)  =  $‘(0;  =  0 


(72a, b) 


J j  (n*®)  ■  •j 


(72c) 


The  solution  is  accomplished  as  before  by  using  a  general  purpose  spline 

collocation  code  to  find  ,  This  is  then  used  as  the  starting  profile  for 
the  unsteady  boundary  layer  solution.  The  details  of  the  procedure  can  be 
found  in  the  Appendix  A. 

4.3  Results  and  Discussion 

As  stated  in  the  Introduction,  the  objective  of  this  investigation  is 
to  compare  the  unsteady  surface  shear  stress  created  on  a  thin  airfoil 
undergoing  transverse  oscillations  to  that  created  by  a  transverse  oscil¬ 
lating  flow  field.  In  both  cases,  a  steady  uniform  flow  approaches  the 
airfoil  from  far  upstream.  Also,  as  discussed  in  the  Introduction,  the 
transverse  oscillation  is  created  by  a  wind  tunnel  with  a  time  varying  wall 
porosity.  Thus, the  actual  spatial  distribution  of  this  velocity  field 
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that  the  airfoil  experiences  will  in  general  be  a  function  of  L/H,M_^,  and 
the  entrance  length  from  the  beginning  of  the  porous  section  to  the  leading 
edge  of  the  airfoil.  Typical  results  of  the  spatial  distribution  calculated 
by  incompressible  potential  theory  are  shown  in  Figure  6.  The  details  of 
this  solution  are  outlined  in  Appendix  B.  In  this  figure,  the  transverse 
velocity  at  the  tunnel  centerline  is  shown  for  various  values  of  L/H,  with¬ 
out  any  obstructions  in  the  test  section.  It  is  obvious  that  the  parameter 
L/H  has  a  strong  influence  on  the  spatial  distribution.  For  L/H  1  the 
distribution  is  approximately  Cos  x  whereas  for  L/H  -►“the  spatial  distribu¬ 
tion  becomes  uniform. 

Subsonic  compressible  flow  distributions  are  shown  in  Figure  7  from  f 

Ref.  [4]  for  purposes  of  comparison.  In  this  case  the  Mach  number  is  an 

I 

additional  parameter,  and  the  two  curves  become  identical  as  -  0  as  ex-  , 

pec ted. 

Once  the  various  parameters,  i.e.  L/H  and  the  entrance  length,  are  fix¬ 
ed  the  spatial  and  temporal  transverse  velocity  at  the  tunnel  centerline  is 
taken  as  an  input  for  the  boundary  layer  calculations.  Obviously,  the 
presence  of  the  airfoil  will  slightly  alter  this  distribution,  however 
in  the  present  analysis,  this  "blockage"  effect  was  neglected.  It  can  be 
argued  that  this  does  not  represent  a  serious  limitation  in  the  present 
analysis  since  other  effects  such  as  tunnel  wall  boundary  layers,  non  uni¬ 
formity  of  the  tunnel  wall  porosity,  three  dimensional  effects,  and  higher 
harmonies  created  by  the  actual  oscillating  flow  would  be  ignored  in  a  po¬ 
tential  flow  analysis  which  included  the  presence  of  the  thin  airfoil.  In 
any  case,  the  present  boundary  layer  analysis  would  still  be  applicable 
since  only  the  free  stream  pressure  gradient  and  edge  velocity  would  be 
altered  and  this  is  treated  as  simply  a  boundary  condition  in  the  present 
study. 
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As  noted  above,  in  the  limit  as  L/H  goes  to  infinity,  the  centerline 
transverse  velocity  becomes  uniform  and  equal  in  magnitude  to  the  wall  velo¬ 
city.  For  this  condition  the  oscillating  flow  field  yields  the  same  un¬ 
steady  pressure  gradient  and  edge  velocity  (and  hence  the  same  shear  stress 
behavior)  as  the  oscillating  airfoil  except  in  the  leading  edge  region. 

This  discrepancy  in  the  L.E.  region  is  due  to  the  mathematical  artifice 
employed  to  remove  the  leading  edge  singularity  and  start  the  boundary  layer 
calculations.  Even  with  the  discrepancy  at  the  leading  edge,  the  relative 
difference  of  integrated  shear  stress  for  the  oscillating  airfoil  vs.  uni¬ 
form  oscillating  transverse  flow  is  less  than  1%.  Consequently,  only  the 
effect  of  a  nonuniform  spatial  distribution  in  the  transverse  velocity  needs 
to  be  considered. 

A  spatial  distribution  given  by  h(X)  =  Cos  X  was  selected  for  compari¬ 
son  with  the  case  of  an  airfoil  oscillating  in  a  uniform  flow.  This  repre¬ 
sents  one  the  worst  possible  conditions  created  by  a  zero  entrance  Length, 

L  ,  and  L/H  near  unity.  In  both  cases  the  frequency  and  amplitude  are 
identical.  Typical  boundary  layer  velocity  profiles  for  the  two  cases  for 
various  values  of  a  are  shown  in  Figures  8  and  9,  at  the  50%  chord  location. 
Notice  that  for  the  oscillating  airfoil,  the  free  stream  disturbance  velocity 
monotonical ly  approaches  zero  as  a  approaches  zero.  This  is  apparent  from 
the  form  of  the  pressure  gradient  given  by  equation  (35)  whereas  for  the 
oscillating  free  stream  this  is  not  true.  Also  notice  as  a  approaches  zero, 
the  velocity  profiles  in  Figure  9  contain  very  little  "overshoot",  i.e.  the 
maximum  value  is  approximately  the  same  as  the  free  stream  value.  This  is  to 
be  expected  since  0  implies  a  quasi-steady  flow  in  the  boundary  layer. 
And,  finally,  notice  that  the  unsteady  boundary  layer  thickness  is  approxi¬ 
mately  equal  to  the  Blasius  boundary  layer,  i.e.  n  s  5.  This  is  a  conse- 
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quence  of  the  scaling  employed  such  that  when  a  =  0[1],  6s/6g  =  0[1]  at  the 
mid  chord  location. 

The  results  of  most  interest  are  shown  in  Figure  10.  This  represents 
the  integrated  non-dimensional  wall  shear  stress  as  a  function  of  the  non- 
dimensional  frequency,  a,  for  both  the  oscillating  airfoil  and  the  oscil¬ 
lating  free  stream.  As  expected  both  curves  go  through  the  orgin,  since 
a  =  0  implies  no  oscillation.  Notice  that  for  large  values  of  a,  both 
curves  asymptotically  approach  a  constant  value  and  the  integrated  shear 
stress  for  the  oscillating  airfoil  is  approximately  20%  less  than  that  for 
an  oscillating  free  stream.  However,  this  limit  must  be  interpreted 
cautiously  since  a  was  restricted  to  be  0[1]  in  the  perturbation  analysis. 
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5.0  APPENDIXES 


Appendix  A 

The  system  of  equations  given  by  (62)  and  (63)  along  with  boundary 
conditions  (64)  and  (65)  represents  a  coupled  set  of  third  order  ndary 
value  problems.  The  solution  of  this  set  yields  the  first  order  ur-^te 
disturbance  for  either  an  oscillating  airfoil  or  a  fixed  airfo’^  v.  .ii 
oscillating  transverse  free  stream.  The  difference  between  the  two 
being  the  unsteady  pressure  gradient  in  (62)  and  (63)  and  the  unsteady  eage 
velocity  in  (64c)  and  (65c).  For  purposes  of  this  discussion,  the 
equations  are  written  as 


111  I  II  I 


H'r  =  ^1(^1.  4'^.  4'^,  -f^) 

(A.l) 

'^i‘  '•'i’  '*'r^ 

(A. 2) 

This  set  of  equations  is  solved  by  the  use  of  a  general  purpose 
boundary  value  code  known  as  COLSYS.  This  code  uses  the  method  of  spline- 
collocation  at  Gaussian  points  and  is  described  in  detail  in  reference 
[lOj.  Only  the  main  "driver"  for  this  code  along  with  the  details  for 
obtaining  the  local  inviscid  unsteady  pressure  gradient  is  discussed  in 
this  section. 

The  four  subroutines  unique  to  the  COLSYS  code  are:  FSUB,  DFSUB, 
GSUB,  and  DGSIIB.  FSUB  is  the  subroutine  for  evaluating  f^  and  f^  at  each 
n  in  the  interval  0  ^  n  rig^.  DFSUB  is  the  subroutine  for  evaluating  the 
Jacobian  of  f.|  and  at  each  n-  And  similarly, GSUB  and  DGSUB  determine 
the  appropriate  boundary  conditions  and  the  Jacobian  respectively. 
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The  subroutine  A  determines  the  coefficients  in  equation  (44) 
which  are  given  by  (40)  and  (41).  And,  finally,  the  subroutine  BC 
determines  the  edge  velocities  and  by  solving  the  perturbation 
form  of  equation  (9),  which  is 


9U, 

-  a _ 1. 

3U, 

+  — ~ 

3-[ 

(A. 3) 


Separating  out  the  periodic  component,  i.e.  e’^,  (A. 3)  can  be 
written  in  integral  form  as 

(A.4) 

This  equation  is  then  integrated  using  the  Gaussian  quadrature  subroutine 
QGIO. 

A  highly  schematized  diagram  of  this  procedure  plus  a  listing  of  the 


driver  program  is  given  next. 


implicit  <^-h,0-2) 

sf-PW,SPi  I  ,.'iFPl ,  xs,  2T,npo  ;R,npn7T 
CnMMON/#?/rtLZ*78 
Cf»!iMr"M/4a/ALPHA  ,  i;k  ,FK  ,  c 
COMilOM/ Ah/AK  ,  A  I 

PpAL***  7€TA(fe)  ,  tRpACF<ai>000)  ,  TUL  (  4 1  ,  7  ( **  >  ,  SF  1  *V  (  4 1  )  ,  fal)  r^F  3^  (ai 

1  )  ,  i»Fl  I  (g1  )  ,.SF?I  (41  )  ,SF3T  fill  )  ,  A4(?1)  ,  A  1  f  ?1  1 
TNTF.CFS  TSPACF  (?8i>0)  ,M(2I  ,  IPA-vT  11  )  ,lTOL  (4) 
fcxTFRNAL  FSli»,OF8llH3,GSy»,oei*OFi,FrT«,FCTT  ,F  k,f  I 
C 
C 

tirnMMa?  .  ... 

1  )=3 
M(?)r3 
C 

ALPrlArl  .? 

FKS.S788 

GKa— 

CALL  A(ALPHA,FK,GK) 
on  5  1=1,21 

b  rtHlTFr6,8)I,AK(n,AI(n 
tt  FORMAT ( 1 X , 14 , ?F 1 4.b1 
xs=?o.n 

DL7  =  .Ti?r>n  . 

C  =  5. 

ZTr.03 

i'45  =  S*1  . 

C 

C 

alffts-o  .ani) 

A«  I(,HT  =  X  s 
IF.TA(  I'tsn.noc 
2FT  A(2)=0.0r0 
ZFT4(x)=o.ono 
ZFTa(4)=0.0l)0 
ZcTAf5)=XS 
Zc  TA (81 =XS 

c 

c  .  . 

IPAn( tl=0 
iPAR(c!)rO 
IPARC^lrO 
IPAR(4lr4 
IPA.^  (51=40000 
lPARfK)=2500 
IPAR (7 ) si 
LL1  10  KsR,ll 
10  lPAtv(Kl=0 

C 

on  20  L=l,4 

LTnL(Ll=L  _ 

20  TOI.  (L1=0.1D-X 


ic 


C 


/<■.  =  ()  .  >11)0 

/.•M  Th  (  5j  )  alpha  ,  AtVIGHT 

iO  F'-.PiVA  T  (  1  HI  ,  '  ALi»HA  =  '  ,F7  .  ,  ‘r  PPL  STf^f  A*1=  •  ,  f  1  0  .  7  1 

;>  S  :V  r  (L,  0 
SSlsO.O 

DO  POO  ISsl.Nb 
^■'JKTTF(b,aO)  7S 

4  0  P-)HWATf///'  -  'iF?.'?//) 

c:ali  cf'LSrs (hCOMP,  V ,  alef  r ,  aptph  i ,  zfta  ,  iPAtv, i  tol  ,  ir i  , Fi  xpf  f 

IFSPACF,  IFl  AG,FSufi,riFSilB,f.SUd,&li3uF,iiOLUTN) 

vJRTTF  (fe,  aF>)DP0  7'7,DP02I 

ab  riiHMAT('  DPUZH  =  ',F10.^,'  DPD/I  ,  F  J  0  .  ?  /  ) 

X=0,0U0 

1)0  100 

CALL  AOPSLW( X , Z,F SPAXF » ISPACtl 
t>F  1  R  (  TD  =7  (  1  ) 

5F?rV  f  IC)  =Z  (?) 
bF  i.R  (  TO  =7.  (  A) 
sF]  l(K)sZ(4) 

SF?T  (IO=Z(b) 

(IC)=ZC(^'> 

x  =  x  +  .  )PH*xb 

100  riiAi  r  f'l/F 
X  =  O.C' 

iK)  11)  rp  =  i,n 

WVIIE(6,  IPO)  X,  SFP?  (  IP)  ,SF3«(  TH)  ,.SF?I  (  I P  )  ,  SF  3  I  ( I  P  ) 
.\iX-».<)PS*xb 
110  C0!'!TIDIIF 

iPU  F')R:lAT(6X,Fb.P«4X, 401^,5) 

X=Xt.OPb*xS 

1):  I  1  XO  1  Ps  I  3  »  ‘4  1  f  P 

1  TF  (H  ,  1  PO  )  X  ,  SFPP  (  IP)  ,  BF  3P  (  IP)  *  :jr  cJ  (  IP  )  f  SF  xn  IP  ) 

x  =  x-f0.f'5*xS 
130  rnMTT'Jilr 

b.skrCbP  +  SF  3P  (  1  ) 

J>5I=S'5  1  +  .SF  31  ( 1  ) 

ZS=Z‘^-*-ri,  7 
<n)u  r.iiijiiDUt 

T  1  F  (iS,  300  )  F5H  ,  SST 

3i)0  pop,.' A  T  (  /  1  X  ,  '  TOT  al  ShFAK  .STPFh.a  =',PF14.4l 

ofC'P  .  ^ 

E''*)' 


,  i  s  A  t:  e , 
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siiP«n  IT  inf  f f  ^ ,  7,  ft 

i-^.FiirTT  kF6l*h  ( 4_H,n-/) 

l'  1  UM  /  Al/SFl«r.SFfH,<bFl],,^FF'l/*^»^^/''F(,?N,i'i^'!l/I 
Ci,;f--N:0:j/A?/OLZ  .  7.'^ 

A  ?/llO  ,  ijf)N  ,  vU  ,  V'J.i 

4  /  AL  PHA  ,  i,K  ,F  K  ,  C 

COA'MGi;/A‘i/i)HFO,i)!Jivi,OPW,J,,3PIO,''*:-'Il,,,HTj- 
CO'-!  Tf'j/ Ab/AR  ,  A  I 

I'lMFNSIOW  Z^b)  tF  (  P)  ,  ?F  1  f,  T  41)  ,F,P?w  (41  )  ,  bF  VHP  T  )  r  SF  (  4i) 
OIMF  jSION  aw f pi )  ,  a  I  (?1  ) 

IF  (ZS.FO.P.OOOGH  TG  1P‘) 

XG  =  y /OSQWT 1 /b) 

II.  =  ) 

1  R  -  P 

ULy^.'^Pbio^xS 
XL  =  1' . 'Mill 
XK  =  yi.-i-Pi  < 

b  IF  (  X  .GF  .  XL  .  AI'UI  .  ■  .LT  .y  RT  GO  TO  10 
xL  =  XGy'n)_x 
X  y  =  X  K  ■*•  '  t  L  y 
1 1.  =  I  L  +  1 

1R=  I  k  ♦  1 

I'll  I  T  G  S 

1  i;  .-.F  1  ,<  x=bF  IP  (  1  L)  +  I  F-F  IP  (  1  k>l  -?F  IK  C  n  )  )  »  (  X-XL  )  /  t  “XL  ) 
GF?Rxs3FPR(  1  H  y(  SFPF  (  1  .-<)  -  SFPR  (ID  )*(  X-XL)  /(  yP-yi,  ) 

SF  1  I  X  =  SF  1  I  ( 1  L  J  +  (  Sr  1  1  f  TF  T  -  SF  1.  T  f  TLM  *  (  X  -  XL)  /  f  y  w-xi  ) 
SFPixSvSF?!  (Tl.)  t  (SF?T  (  i  P)  -SFPT  (Tl  )  T  *  (  X-Xu)  Z  (  yp->l  ) 


lFfyF-s.0)4(',ar;,S(j 

4<(  F  •/—  ,0  0b<?PS*y-F  +  (i,  1  ub7P'4*xb*3yp  +  0,vll  44  *  X  P  »  *  X  ••  0 , 0  Ob^)  C:u*XF<t*  —  + 

l<.)*Xh**5 

F  P  =  0  ,  ■<  P  4  1  7  7  *  X  F  +  1) ,  1  b  o  'X  P  *  X  s  *  *  ?  -  0  .  M'.  c,  I  ^  I  «  X  R  *  *  X  -  r  . ' '  0  (■  F  s  4  *  X  '-  *  »  c 
ll4p*Xf<**S 

Fi-'P3(T,'J5Pub”'0..  014if47*X;'t  +  0,0Xb4^'i*X4**P'-ii,(iX7-,iio*x^*»7*r'^:,; 
l*4-0.0  0  0bXbA*xF>**'> 

f.n  Tii-bO 

bO  F0=X.PF3Po+(  Xr'-S.o) 

KP= 1  .  0 
Fop  =  n.r> 
h  ('■  11  n  3  F  P 

iji,iN  =  Fpr/2b**n.b 

V '131. 0/73**  !I.‘=*(XO*FP-FO) 

V'G'‘  =  0.5/7S*yB*FpP 

c 

c 

If  (/F.to.  .OPOOTGO  TO  XOi) 

T  NHk  =  ',I0*3FPw  X  /hL/-'-)MA'*SP  1  WX  /:'L7 
1  ,v)NI=UU*.*lF?T  X/DL7-!tGiy*SFl  TX/.Tl./ 

zc=l.>i/(?.n*Z.s**t.5*o.sr7FUi.o-'^')) 

i'PPZW  =  7r  *AW  (  1  ) 
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Slipi^'n  ITIMF  F (  X  ,  ?  ,  F1 
I'-.PiirTT  FF  Al.  *><  (  A -H,  n-/ ) 

H;f--ivo;j/A;5/0LZ  ,7.'^ 

Cn^■.v;n^>/ A  ?/iiO  ,  ijOr  ,  vu,  ViJ.-j 
Cni-MO'-/ Aa/  ALPHA  ,  I,A  ,Fi<  ,  C 

/  ,npp  1 ,  OPw;j ,  ,)P  1 0 ,  ‘'--M  1 , 1'i-’i  ? 

CO'-!  .IP'J  /  Ab/ AS  ,  A  1 

‘MMFAJSTOM  Z(b)TF  (  P)  ,  SF  1  f.  f  i*n  r  (^n  ,  U  f  p  n  r  SF  PI  (  Ui) 

OIMF  jSIOW  aw  f  PI)  ,  A  1  (?n 
IF(7.s.FO.O.O()0)G(«  T-0  IP-) 

X =  V  / 0  s Q PI  I  / S  ) 


II.  =  ) 

)  ■•'  -  P 

ULX  =  .0PS'^0*»S 

XL  =0  .  'MM’S 

X  K  =  y  (.  .  ■•■  P I  X 

I  F  f  X  .';F  .  XL  .  AOu  .  «  .L  T  .  V  w)  GO  TO  I'l 
xl  =  Xl+‘)LX 

X  .-V  =  X  K  •*• '  I L  X 

I I  =  I  L  +  1 

IPrlw+l 

III  I  1  '!  S 

1 1,  .'iF  1  X  =bF  T  p  (  T  L  )  +  I  ^,F  IP  (  1  P  ) -?F  IH  C  Tl  )  )  *  (  x  -  Xu  )  /  (  VF  -  x  l  ) 
bFp.y  x=3Fp-y  ( 1  L  )  ♦  (  SFPF  (  1  O  -SFPR  (  ID  )  *  (X-^L)  /  (  XF-XL  ) 
.Sr-  1  I  X  =  SF  1  I  (  1  L  J  +  I  3r  1  U  T?  1  -  SD  I  f  Tin  *  (  X-  Xl.)  /  f  >  S'-Xl  ) 
SFPIxsSFPI  ( It, )  ♦  (  3FpT  (  W)  -SFPT  (  Tl  ) )  *  (  x-Xl)  /  (  XK-XL  ) 


IF  (  XF-S.  n  )  6c  ,  ao  ,  S(j 

Ul;  F.ji.0()FP?S*XH-K).TUb7qO*Xb*ytP  +  0.«>’H13^*)‘P*»?«'f>  .  0  f.'buJ  (i  u  *Xli>  *  u  t  i;..  u  ti  1.' U  1 
lF*Xh**5 

f  P  =  0  ,  P  a  1  7  7  *  X  F  0  ,  (/  1  b  <'  3  P  *  X  P  *  *  ?  -  0  ,  '  I'l  V  i  -S  1  *  X  p  *  *  X  -  f;  ,  '  I  F,  S  9  *  X  »  u  -f  i';  ,  ( I  I'l  ( I  r 

}ii*P*Xii**s  - 

Fi-'i^  =  ii,3X?ub”0.O1Fpii7*XO+’f',OXbt;f''i*X’''**p”'i,0^7^''9'AXF. 

1  HI,  -  0.0fi0i'3b‘^yrXF**s 

r.n  T;j-bO  -  . 

bO  Fi)=  X.PriPo+ (  xr-s.n) 

FP=1  ,  rt 

Fopzo.n  .... 

bi'i  iiO  =  FP 

Ui,rM  =  Fur/;jS**n  .F 

VO=l,')XP,G/7ii**(i,p*fXP*FP~F0)  — 

VOi''  =  0.5/7S*Xrt*FPP 


If- ( /s.e  0 . . upfioi  GO  TO  .xoo 
T  i.Hfyr'iOiySFPS'  X  /OL/-'.J|)^'*SF  1  RX/i>L7 
T  ,v!HT=|.iO*RF?T  X  /DL7-m>iX*SF  \  IX  /.)L7 

zc  =  i.  *)/(?. 0*7. s**i. 5*  DsrjFTd.o-'^')) 

iinr.7R  =  7r*AR  f  1  ) 
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r.it>p/]  =/c'  *A I  ( 1 ) 

UO  70  Ms?, 10 

XMSV.-1 

A7  =  XM*»)«R|[;i)Sf?.(i*7  5'-1  .0) 

JiCsa  .  0»  ?S*XN;*Or()S  (  A/ ) 
f .  P  n  7  R  =  n  P  D  7  K  -  S  C  *  /'  r  *  A  K  (  M 1 
7t.  I^pn  /  T  =(»Pn7  T  -  SC  *7C  *  A  I  fr’l 

OFIPsDPDZI 

(■  (  \  )=vn*Z  C^)  -  (VC'^'-Un/OL?)  *2  <  1  )  /OL7.-ALP^A*2  f  Sj  ♦-.POZ-f  -  TiHfV 

Ffpi=vu*2(M-(vrif,-iin/aLZi*7(S)-L..'->7(0/nL7tAi.^HA.7,?)..rP, 

R  t  T  I J  R I ) 


IPO 


pno 


Piu 


a?o 


X^^.rP.*X*nsjRT  (Cl 
IF(xb-2.0)  ?n(>,PO()r^l  0 

F0  =  -.  001  07*XH  +  .(rPa32*XP*X  .  I  7  *  XP  n- *  7  . 
FP=l!?<t*XR  -  .SX*yp*XH  +  .0F4*XH**X. 

c;;i  rn  ??o 

Fn=i.ih?a  +  .ppp*(xR-p.O) 

FP=1  .0 

\;Os-?.OoFO*nSUWT(C) 

Tia' r -  'i  ,  0  *F  P  »  C 
l.)Pt  7F  =  - alpha  *  ALPHA  /  ?.  0 
L)Pn/ 1  s?  ,0*  Al.  PriA*C 
OPR.jsliPLZH 
iiPT0  =  |)P07I 


30y 


31  R 


♦  .i.,?pac,*xB**4, 


F  f  t  Isi/O *2  (  3)  -  7  (P)  -  A I  pH A*7,  (5)  +'^P''  ZP 

F(P)sVO*ZtS)-VOF;*2(‘')+AlPHA*7(P)+‘'P'7I 
HF  I  u«rj 


ZC.  =  1  .0/  (P.O*Z  (*-*  '  .i*l).SUPT  (  1 ,0-71 
HPOZP-ZC* AP ( t ) 

OPDZ I  =  ZC  *  A  T  f 1 ) 
l>i  I  X  \  U  N>.sp  r  1  0 
X  M  =  -  1 

A  7  =  X  M  *  n  A  !  ?  C  O'  3  f  P  .  0  *  /  T  -  1  .01 
bC=a.0*7T*XM*t;(;uS(  AZl 
0  P  0  7  R  =  0  P  0  7  F  -  S  C  *  7  C  *  A  «  (  A. ) 
OFOZT=l)POZT-vSC*vr*AI  IMl 
.UPDZT=(nPUZT4UPI0) /2.0 
UPD7P=(npr)ZR+OPP01/2.0 
pop  1 sDPPZR 
uPTisOPnzi 

VrviHR  =  iin*SF2FX/rt|.  7-il(iiV*SPl  RX/nL7 

T>JHTr'jO*SF  ?1  X/Pt  Z-UOM*SF  1  1  X/I)L7 


F-(1  jrvn.Zf  ;l-(VO^-'Jl.'/OLZ)*Z(P)-''n^'*Z(l)/0L.7-ALPF'A*7(^,)+^P■'7H-^^^'^ 

•  iUiM*7  (  «) /OL  Z  + AI.Ph  A*  Z  (  2)  ♦OP ''Z  1-T  Ni-il 


F  (P)  =  Vf'»Z  (6)  -  I  Vl'f  -tjO/OLZ)  *7  ' 

RFTUW'J 


SIlh-^OliT  OFSmWV, /T.iif  ) 

l^'F'LirTT  RE«L*f<  (  A-H,ri-/) 
COMMO.V/ A3/F>LZ  »  7  8 
COMMON/ A  ?/Ur; ,  L'ni,i,  VO, 

CnMMON/A«/ALl^HA,GK,FK,f: 

0  I>‘f  Nil  (*A:  0F/P',<'),7(A>) 

c 

IF  (7  8.^0.0.0  100  TO  ?« 

0F(  1 , I )=-0DN/uL7 
OF  ( 1  ,  <?T  =-vON  +  Of>/l<L  / 
i;F  (  1  ,  =Vfi 
-uM  1  ,<*)  =0.000 
UF  (  T ,5) =-ALPHA 
UF  f1  ,N) =0.000 
OF(i»,  \  )=0.«00 
OF  (,?,  P)  =^0?^^  A 
0F(?, ?) =0.000 
(',p  (2 , 0  )  =-UO-n(/UL  / 

[)(■'  f  ?,S)=-VOiM  +  illT/t)L/ 

OF (P, 01 = VO 
►  T  1 1 W  N 
C 

?0  OFd  ,  11=0.000 
JF  (1  ,P)=-V.)N 
OF  f 1 ,  T) =V0 
OF (1 ,4) =0.D0 
OF ( 1 , S1=-ALFHA 
OF  (  1  ,ft)  =0.0110 

nF(p,n=o.ono 
i)F(p, p;=ah^na.  . 

OF(P,  ilrO.OOO 
OF (P, 41 =0.000 
i)F  (  p ,  ‘>1  =  - vnN 
OF  CP,h)=VO 

K  F  T 1 1  k  fs( 

END 


SMFKO  iTINr:  •- SOd  (  T  ,  7  ,  Gl 
TMPLIFIT  :<F  a  L*  B  i  A-H  ,  0-7) 
);[Mf..vSinN  Z(fa) 

c 

r.n  TO  (i,?,7.,4rB,f'),T 
1  ti=2iu 
kETUPV 
P  G=Z(?1 
KFTUB'.V 
3  G=7(4) 

I^FTUR!'! 
tt  G=Z(5) 

PFTI'R'J 

d  CALL  AT  (D  -e,  0)  1 

0=2  (2) -'Jp . 

KFTMP 
h  CALL 

.i  =  Z(S)-Ol 
kf-  T  I'k'. 
e  '\ID 
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S'JPkUllT]  I\li  nGSUM  (  I  ,  Z,OG) 
implicit  f?FAL  •'*  f  A-H,0-Z) 
UiMEN-SIOi'i  2fb3,aG(ft). 
un  10  J  =  1 » 6 
10  GG(J1  =  O.OPO 

.  GO  TO  (1, 2, 0,5, 6) ,  I 
1  iJG{n  =  i.coo 
WFTiJP'J 

i  DG.C2)  =  1.000 
KKTuOi'J 

3  UGf*4)  =  l  .01)0 

M  E I U  w  IV' 

0  aG(5)  =  l  .Oix) 
h<FTl|PN 

5  OG(2)=l.OnO  . 
ketiipm 

fc  0G(5)  =  1  .00(1 
PETilPN 
feivin 


oUBkODTINF  ft f AL Hh a , F K , GK ) 

IMPLICIT  PEfL*a(A-H,G-7) 

DIMCrijIUN  AH  C21I  .  Al.L^n  ,PkC2?)  .  °I  (22) 
CnN‘,,\iON/Ah/AH,  AI 
COM, 1PM/ A  7 /K  1 
fcXTEPNAL  FR,FI. 

00  50  N=l,?2 
K  l=N' 

XL=U.Q  . 

Kti=?. 1415927 

CALI.  OG10(XL,X(j,FR,PMP) 

CALL  GG10.(XL,XU,FI,PMT) 

PR  (,vi)  s-PNP/3. 1  41  5927 
50  Plf''i)=-PNT/3. 141  5927 

ARf  n=FK*,IPR(  J).+  P9(2)  )-PR(2)  -  GH*  ( .2  1  f  U.+F  I  (21  )  .. 
A1  (1  )=(;k*  (HR  (1  )+PK  (?)  )  +  FK*(PI(1  1+P]  (?))-H\  r?) 
(>0  60  M  =  2,21 

AP(N)_=-.1.25*ALPH4/CA'-1)  )*..fPI  (N-U-PJ  (1^+ l )  } +H2 (Ji) 

60  AT  (N)  =  (  .25*  alpha  /  (N-n  1  *  (PK  (N-1  )  -PR  (^+1  1  )+P  I  (f  ) 

RETllRfM 

two  .. 
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sut^-yOijr !fJt  i=<C(Ut?,un 
Ii'^PLTriT  RFAL*P  (A-H,0-7) 

Cn^MOM/AP/Ol  Z»7S  . 

CnMMOiJ/A4/ alpha  ,RK  ,FK  ,C 
CnviflON/AS/OPHOiDPRtl  ,DPI? 

tXTePN-AL  PCT^-,F€TI 

lF(7.S.Lrt.0.(''5l'0)Gl»  TO  10 

XL  =  f).04 

XU=7S 

CALL  e&10(-J^.,)ty-rFCTfc^,VK)~ 

CALL  fjG10(XL,XL!,PCTI,  VT) 

UR  =-Yf<*OCOS  {  AL  PHA*7S.l  -  Y  I  *0S  (  U,  Pni  »2  ^  UR? 

tLi  =  yr*omvlalpm*zs)--  Yi*oc:as(ALH’uA<-z.s)  ♦  -u-i^ 
RFTijRf- 

JO  If  (ZS-0.0pni;)a4>,io,?<' 

?.{)  URs'JRJ  -!JL7*  (uPRl+i-'PR?)  /P.O 
UI=(JI1-DL7*  (UPIl+nPIPl/P.O 

UC?  =  UR  - .  -  . 

UIP  =  UI 
Kt.  TURN 

XO  UPaURO-Ol.Z*  (uPf'<U»i)PRt) /?.0 
iJIsuin-OL/*  fl'PIo+UPll ) /^.n 
UPI=I.JR 
I'llsiJT 
HeTliRN 
00  UP  =  ().0 

U^s-A-LPHA/^.O  .  ... 

UROsUP 

UtOs'^I 

*<FT(»Ri'( 

ewn 


.‘illHRCiilT  I  JF  fjf,  1  0  (  XL  ,  XU,FCT,  V  ) 

ImPlTCTT  RFAL*H  (A-m,(!-7) 

EXIEWWAL  FX^TR^FCll  ,FR,FI 

AsO.S* (XutXL  1 

osXU-YL 

C=0.40fe9535*rt  -  - 

YsO.O  4X<X5b7* (FCT ( A+C)+FCT ( A-r1 1 

C  =  0.4X?‘S3J  7*0 

YsY.*.07«7i567*(f  CTiA+n.)  ♦FCHA.-Cl) 

C  =  0.3  AP7f)«ft>*8 

YsY  +  O.  }  00543?*  fFCT  fA*C)-»-FCT(A-r)J 

CsO,aiF.<»077*a  . 

YsY  +  0.1  346334*  (FCT  (A ♦D+FCT  (A-D  ) 
C=0.07a45717*b 

Yse*  CY  +  0.1  '477fe2l*fFCT  (A*C)+FCTCA-C)  1  ) 
RETURN 


I 


K'lr. CT  rr.r>!  i^kT’IF  IA) 

rM->L  rri  I  AL*f'(  A-ri,(i-/) 
COM;40w/Aa/ALPHA,G)'  ,FK,C 
C.PM''irn-i/ A7  /K  1 
VT  =  ALHhA/2.(t 

F  rrVI*OCOS(  (K  1-n  *TnFTA) 

P  F  T 

EMO 


Fi'’'!CT  [t>i,  j.-prn.Fi/i) 
iMPl.TCTT  AL*^  (  A-H,  (i-7  ) 
CnMMPN/A7/K J 
VP  =  0.0 

Fp=\/k'»nciis  (  f  K  1  -1  )  *,TH,F  T  A  ) 
kF  1  ilkM 
C  W  0 


r  MHCTTni','  FCTTfZS) 

I  f-'PL  T  C  T  'I  kt*  A  L  *  P  f  A  -H  ,  0-  7  ) 

OIHEiVSION  Ak.i  21)  ,  Ai  (?1  ). 

Lf'^  >«/ ALPHA  ,  GK  ,FK  ,  C 

C('m.«/iOM/aF/AR,  A  I 

/C  =  l  .0/(2.0aZS**]  .5*0SaRT(Ui.*-2ij) )  .  .. 

FrTT=f'STW(ALPhA*ZS)^7C*AR(n+r'(.'-S(ALPrtA*7S)*7.i:AAI(n 
•)n  tO(; 

X’-'sH-l  , 

tt7sxvt*r'Akrns(?.n*7s-i  .o) 

.  I)*  7S«  (  A  7 1 

1  Du  t-CTIsFr  T  T-OS  I  A.  (AlPMA^Z*))  *ZC*bL*Ak  f  ">  J  -OC  D  S  I  A  lF  k  A  »  Z  ^  ZC »» i  r  1-  ) 
r<  F  T  I  Ik  '  'I 

c  T' 


FilAriTflA'  FCIkfZS) 

l-ipi.jrTT  kFAL*«  (  A-H,(,-7) 
tJTMFN?,irjN  Aw(2n,Ar(?l1 

C0i«’Mn^/Aa/ Alpha, GH  ,fk  ,c 

CuM  -^OiJ/Ah/AK  ,  M 

Zf;5:l.i>/(?.0*73**l.'5*USfjRT(1.O-7.s)) 

»-CT'V  =  i'>Ci'S(ALPHA*7?)*7C*AP(n-i>S['>i|A|.PHA*Z;?)*7r*AT,  rn 

T'l  IDO  ''  =  2,10 

XMrM-1 

A7  =  XNi*r-APC0S  f2.o*7?*-l  .0) 

.sr  =  D.O*7S*XM*L‘CO.S(A7) 

100  PCI  ksFCTa-OCOF-f  Al,FkA*7S)  *Zr.*.Sr*AP  trrt)  +nsi  N  (  ALPHA*/h)j>ZC*SC*AI  fM 
kE  TijtJfj 
t'ln 
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The  actual  flow  field  induced  by  temporal  variations  of  the  test 
section  wall  porosity  to  an  otherwise  imcompressible  uniform  wind  tunnel 
flow  is  a  function  of  the  Reynolds  number,  entrance  length  (Lg),  L/H,  and 
the  geometry  of  the  model  in  the  test  section.  Other  effects  such  as  free 
stream  turbulence,  the  details  of  the  flow  through  the  porous  section  and 
the  frequency  of  oscillation  can  be  important.  However,  a  reasonable  ap¬ 
proximation  to  the  unsteady  oscillations  can  be  obtained  in  analytical  form 
if  a  potential  flow  model  is  adopted.  This  is  the  same  basic  restriction 
employed  in  reference  [4]. 

Refering  to  Figure  11,  and  considering  the  effect  of  the  variable  wall 
porosity  without  any  obstructions  in  the  test  section,  the  potential  function 
is  given  by  solving 

V^(j>  =  0  (B.l) 

subject  to  the  following  boundary  conditions 

Vg(x)  Cos  “t  ,  -L<x<l^ 

0  ,  1]  >1  (B-2) 

z 

->  and  1^  -*0  as  jx|  -*■  «>  (B.3) 

The  solution  to  this  system  of  equations  can  be  found  in  closed  form 
by  using  the  Fourier  Transform  Method.  First  writing  <J)(x,y,t)  =  (J)5(x,y)  Cos  w 
as  a  consequence  of  linearity  of  the  system,  and  then  introducing  the  Fourier 


1 


transform  defined  below. 


4>(s,y)  =  41  (x,y)e"’  ^  dx 

/  2v 


The  original  system  of  equations  becomes 


(B.4) 


d^<I>  A  -  n 

— 2  -  S  <!>  -  0  , 

dy^ 


(B.5) 


with  the  boundary  condition 


'y=0,H 


V^(s) 

o'  '  •> 


(B.6) 


where 


Vo(s)  =  ~  V^{x)  e*'^^  dx 
’^2tt 


Solving  for  d4>/dy,  which  is  the  variable  of  interest,  we  obtain 


^  _  y  /  X  Sinh  sy  +  Sinh  s(H-y) 
dy  *o'^'  Sinh  sH 

Taking  the  inverse  transform  and  using  (B.7)  yields 


(B.7) 


(B.8) 


+  1/2  f+< 


04',  1 

IT'S  Vt’ 


Sinh  sy  +  Sinh  s(H-y)  is(x-5) 
Sinh  sH  ® 


ds^  d^ 


(B.9) 


Integrating  (B.9)  with  respect  to  the  transform  variable  s,  yields 


where 

A"  =  Cos  ^  Cosh  ^(x  +  ^) 

B-  =  Cosh  ^(x  -\)  t  Cos  5^ 
C*  =  Cosh  J-(x  +  |)  1  Cos 
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Figure  2. a.  Oscillating  airfoil  in  a  uniform  stream 


I  VQ=-awSinuit 


Figure  2.b.  Fixed  airfoil  with  an  oscillating  transverse  flow 


Figure  3.  Potential  flow  and  boundary  layer  coordinates 
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Figure  4.  Conformal  mapping  of  line  segment  into  a  circle 
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Longitudinal  distribution  of  flow  angle  amplitude  for 
a  symmetric  ball-wall  installation,  L^H  =  1 .0. 
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Longitudinal  distribution  of  flow  angle  amplitude 
for  a  symmetric  ball-wall  installation,  L/H  =  3.0. 


Figure  7.  Spatial  distribution  for  subsonic  compressible  flow  from  reference  [4] 
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Figure  10.  Nondimensional  integrated  shear  stress  vs.  nondimensional 
frequency 
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Figure  11.  Geometry  and  notation  for  the  potential  flow 
solution 


